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ABSTRACT 

The observational signatures of the first cosmic explosions and their chemical imprint on second- 
generation stars both crucially depend on how heavy elements mix within the star at the earliest 
stages of the blast. We present numerical simulations of the early evolution of Population III pair- 
instability supernovae with the new adaptive mesh refinement code CASTRO. In stark contrast to 15 - 
40 Mq core-collapse primordial supernovae, we find no mixing in most 150 - 250 M pair-instability 
supernovae out to times well after breakout from the surface of the star. This may be the key to 
determining the mass of the progenitor of a primeval supernova, because vigorous mixing will cause 
emission lines from heavy metals such as Fe and Ni to appear much sooner in the light curves of core- 
collapse supernovae than in those of pair-instability explosions. Our results also imply that unlike 
low-mass Pop III supernovae, whose collective metal yields can be directly compared to the chemical 
abundances of extremely metal-poor stars, further detailed numerical simulations will be required to 
determine the nucleosynthetic imprint of very massive Pop III stars on their direct descendants. 

Subject headings: 



1. INTRODUCTION 



The first stars in the 
are l ik ely very m a ssive 
[1991 lAbel et al.l I2000L 



universe form at z ~ 20 and 
30 - 500 M(7) (Bro mmet all 



20021 iBromm et al.l |2002l 



Nakamura fc Umemural 120011 : lO'Shea fc Normanl l2007h 
The fates of these stars depend on their masses: 15 
- 50 M stars die in core collapse supernovae (CC 
SNe), 140 - 260 M stars explode in far more ener- 
getic thermonuclear pa ir-instability supernovae (PISNe 
iHeger fe Wooslevl l2002f ). and 40 - 60 M stars may 
die as hyper novae, whose explosion mechanisms are 
not yet understood but are thought to have energies 
i ntermediate to thos e of CC and pair-i nstability SNe 
(llwamoto et aLll2QQ5h- Most Pop III SNe (IBromm et al.l 
120031 : iKitavama fe Yoshidal 120051 : iGreif et al.l l2Q07h oc- 
cur in low densities (0.1 - 1 cm -1 ) because UV radia- 
tion from the star sweeps most of t he baryons from the 
dark matter halo in wh i ch it resides (IWhalen et al.ll2004l: 
I Kitavama et al.l l2004t lAlvarez et al.l 120061 : lAbel et al.l 
120071 : IWise fe AbelH2008aT ~ 

Metals from Pop III SNe determine the character of 
second-generation stars and the primeval galaxies they 
populate by enhancing cooling in the gas in which such 
stars form, and therefore the mass scales on which it frag- 
ments. The manner in which the first metals contami- 
nate pristine gas in the primordial IGM crucially depends 
on mixing processes that have only begun to be studied 
numerically. Preliminary calculations indicate that met- 
als from Pop III SNe mix with gas in a halo on two 
characteristic spatial scales prior t o their emergence int o 
cosmological flows on kpc scales ([Whalen et al.l l2008[ ): 
10 - 15 pc, when a reverse shock forms in the remnant, 
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and 150 - 200 pc, when the remnant collides with the 
dense shell of the relic H II region of the progenitor. 
New simulations now prove that violent mixing can occur 
within the star prior to s hock breakout from it s surface, 
if the star is 15 - 40 M (|Joggerst et al.ll2Q10aL hereafter 
JET10). Capturing mixing on the smallest scales is pre- 
requisite to following its cascade out to the largest ones, 
and hence to determining the colors and morphologies 
of primitive galaxies. Early mixing also governs which 
elements are imprinted on next-generation stars, whose 
chemical abundances may impose indirect constraints on 
the masses of the stars that enriched them. Low-mass 
remnants of this generation are now sought in ongoing 
surve ys of ancient, dim m etal- poor stars in the Galactic 
halo ([Beers fc Christliebl 120051 : iFrebel et al.ll2005h . Since 
mixing in Pop III SNe in part sets elemental abundances 
in second- generation stars, it is integral to future mea- 
sures of the primordial IMF. 

Early mixing is also key to computing the light curves 
and spectra of the first cosmic explosions, whose detec- 
tion could yield the first direct measure of the primordial 
IMF. Optical and UV radiation breaks free of the SN 
shock when the ejecta reaches the outer envelope of the 
star and is exposed to the IGM. In the frame of the shock, 
the photosphere from which photons escape into space 
descends deeper into the ejecta as it expands outward 
because of the spherical dilution of the ejecta. If mixing 
precedes radiation breakout, it determines the elements 
that the photosphere encounters as it sinks deeper into 
the ejecta, and therefore the emission lines that propa- 
gate into the IGM over time. Accurate multidimensional 
models of the explosion from its earliest stages are there- 
fore necessary to compute lines in primordial SN light 
curves and spectra. 

Mixing in galactic core-collapse supernovae has been 
studied for over twent y years, particular ly in connection 
with SN 1987A (e.g. iFrver et al.l I2007L and references 
therein). A core objective of these studies is to under- 
stand the premature appearance of 56 Ni in the spectra 
of 1987 A, whose emission lines appear much earlier than 
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a simple picture of segregated, spherically-symmetric ex- 
panding mass shells would predict. Mix ing in Pop III 
core-c ollapse SNe was first examined by iJoggerst et al.l 
(|2QQ9h . who found that both vigorous mixing and fall- 
back onto the compact remnant in 15 - 40 M Pop III 
SNe govern which metals escape into the IGM at high 
redshifts. 

PISNe may have been commonplace in the primeval 
universe, but their enrichment of the early IGM is yet to 
be understood. To investigate the propagation of met- 
als into the IGM by such explosions, constrain their nu- 
cleosynthetic imprint on second-generation stars, and to 
evaluate the impact of mixing on PISN spectra, we have 
performed two-dimensional simulations of the explosions 
of 150 - 250 M stars. In § 2 we review how blast profiles 
from the KEPLER code were ported as initial conditions 
to CASTRO and discuss what factors govern the presuper- 
nova structure of the star. We examine the outcomes of 
the explosions in § 3 and in § 4 we conclude. 

2. MODELS 

2.1. KEPLER 

As in JET10, the simulations in our PISN survey were 
carried out in two stages. First, primordial stars were 
evolved through all stages of stable nuclear burning from 
the zero-age main sequence to initial collapse via the pair 
instability in the one- dimensional Lagran g ian stellar evo- 
lution code KEPLER (|Weaver et al.l 119781 : IWooslev et all 
l2QQ2h . The PISN is triggered when this collapse incites 
explosive O, and some cases Si, burning. Unlike the 
simulations of JET10, in which the blast is artificially 
launched with a piston, the pair instability and subse- 
quent collapse that triggers these explosions are an emer- 
gent feature of the stellar evolution calculation. They are 
genuinely spherical, barring (magneto) rotational effects, 
and their energies are set by O and Si burning. The blast 
was followed until the end of all nuclear burning, ~ 20 
s after the start of the explosion. The energy generated 
was computed with a 19-isotope network up to the point 
of oxygen depletion in the core of the star and with a 
128-isotope quasi-equilibrium network thereafter. 

2.2. CASTRO 

The one-dimensional explosion profiles were then 
mapped onto two-dimensional RZ axisymmetric grids 
in CASTRO (Compressible ASTROphysics), a multi- 
dimensional Eulerian AMR code with a high- ord er un- 
split Godunov hydro solver (Al mgren et al.ll2010l ). Each 
explosion was evolved past breakout from the surface of 
the star until all mixing ceased and each element in the 
ejecta was expanding homologously in mass coordinate. 
We smoothly join the density at the surface of the star 
to a uniform circumstellar medium of 1 cm -3 with an 
r -31 power law, in keeping with the usual assumption of 
a low-density H II region around the progenitor star with 
no wind-blown shell or prior mass ejection. The medium 
beyond the star has no effect on the dynamics within the 
star if its density falls more steeply than r -3 . In mapping 
the radial profile onto the RZ grid in CASTRO, care was 
taken to resolve the key elements of the explosion: the 
shock, the shells of elements, and the high-density core. 
In particular, both the 56 Ni core and the O shell were re- 
solved with a minimum of 16 cells. Our mapping excludes 



departures from spherical symmetry due to O burning, 
but such perturbations would likely be high mode and 
low amplitude and therefore have minimal effect on the 
evolution of instabilities. Our models thus only capture 
later asymmetries of mode greater than / = 1 or 2. 

As in JET10, we adopt a monopole approximation for 
self-gravity. We first compute a radial average of the den- 
sity from the RZ grid to create a one-dimensional den- 
sity profile. We then compute a one-dimensional gravi- 
tational potential from this profile and map it back onto 
the RZ grid. Since departures from spherical symme- 
try in the densities are minor, this approximation has a 
negligible effect on the final state of the explosion. Be- 
cause the PISN completely disperses the star, there is no 
compact remnant, fallback, or thus any need to include 
a point potential centered at the origin of the coordinate 
mesh. 

We follow 15 chemical elements as individual species, 
each with their own continuity equation, and calculate 
local energy deposition due to radioactive decay of 56 Ni 
in the same manner as JET10. However, the formation 
of a nearly degenerate core at white dwarf densities in the 
progenitor necessitates the use of the Helmholtz equation 
of state (EOS) at early stages of the explosion. As the 
ejecta expands and cools, we transition back to the ideal 
EOS used in JET10, which assumes that the gas is fully 
ionized and includes contributions from both radiation 
and ideal gas pressure. 

The base grid is 1024 2 , with the initial outer bound- 
aries set so that the inner portion of the star is resolved 
as described above using no more than 6 levels of refine- 
ment. The star is centered at the lower left corner of the 
mesh. We apply reflecting and outflow boundary condi- 
tions to the inner and outer boundaries of the grid, re- 
spectively. Our refinement criteria are the same as those 
in JET10. When the shock nears the edge of the grid, 
the simulation is stopped, the grid is doubled, and the 
calculation is then restarted, subtracting or adding lev- 
els of refinement as needed. This procedure is repeated 
up to 12 times, depending on the model. We halt the 
simulation when all chemical species are expanding ho- 
mologously, which always occurs by the time the ejecta 
has propagated a short distance into the uniform circum- 
stellar density. 

2.3. Progenitor Models 

JET 10 found that mixing in low-mass Pop III core- 
collapse explosions primarily depends on the internal 
structure of the progenitor. We likewise expect the early 
evolution of pair-production explosions to be determined 
by the envelope of the star, which is determined its mass, 
internal convective mixing over its lifetime, and by its 
metallicity. Capturing the full range of structures for 
these stars is essential to a comprehensive survey of early 
mixing in Pop III super novae. 

2.3.1. Convective Mixing 

The initial absence of metals and the large contribu- 
tion of radiation to the pressure in massive Pop III stars 
promotes convection within them. The CNO cycle can- 
not begin in primordial stars until a threshold mass frac- 
tion of 12 C is first created by the triple a process. This 
trace 12 C sets the entropy of the hydrogen layer to be 
just above that in the core, without the sharp entropy 
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TABLE 1 

PISNe Models: Properties at time models were mapped to 2D 



model 


7\/f f~\K \ 


M N (M©) 


M56 Ni (Mq) 


M final ( M 0j 


ri (1U CII1 J 


&kin (10 ergsj 


ul50 


41 


9.1(-5) 


0.079 


143 


16 


3.7 


ul75 


46 


l.l(-4) 


0.72 


164 


18 


9.5 


u200 


50 


1.3(-4) 


5.2 


183 


19 


17 


u225 


54 


l.l(-3) 


17 


200 


33 


28 


u250 


68 


1.7(-4) 


38 


238 


23 


34 


zl75 


53 


1.6(-5) 


0.24 


175 


4.2 


4.7 


z200 


60 


1.8(-5) 


2.0 


200 


4.5 


12 


z225 


67 


1.7(-5) 


8.8 


225 


4.9 


17 


z250 


74 


1.6(-5) 


23 


250 


6.2 


29 



gradient in the upper layer of the helium shell that is 
usually present in He burning stars with metals. This 
plus radiation pressure facilitates convection. In 140 - 
260 M stars, the central convection zone can approach, 
come in contact with, or even reach into the lower hy- 
drogen layers, mixing them with carbon brought up from 
the core during helium burning. When these two high- 
temperature components mix, they burn vigorously, ele- 
vating energy release rates in the H shell by up to several 
orders of magnitude (the so-called hydrogen boost). 

Convection affects the structure of the star in two ways. 
First, since it raises energy production (and hence opaci- 
ties) in the lower hydrogen layers, the star can puff up by 
more than an order of magnitude in radius and acquire 
a red supergiant structure. Second, if convection is ex- 
treme the transport of material out of the core could re- 
duce its size and explosion energy in comparison to mod- 
est convection. Unfortunately, one-dimensional stellar 
evolution models cannot predict these inherently three- 
dimensional processes from first principles. Instead, they 
parametrize them with semi-convection coefficients. 

2.3.2. Metallicity 

Gas in high-redshift halos that is enriched to metal- 
licities below 10 -3 ' 5 Z fragments on mass scales that 
are essentially identical to those of pristine gas and 
still forms very massive stars (e. g. iBromm "eTa!][200TI : 
iMackev et all 120031 : ISmith fc Sigurdssonl 1200/D . How- 
ever, such small metal fractions are more than enough 
to enhance CNO reaction and energy generation rates 
in the hydrogen burning layers of the sta r, enlarging it 
in the same manner as convective mixing (jHirschil l2QQ7t 
lEkstrom et al.l [2QQ8h . Hence, we would expect a strong 
degeneracy between the influence of metals and convec- 
tion on the envelope of the star, and that the full range of 
mixing in PISNe can be as easily spanned by metallicity 
as by convective overshoot. 

We considered 150, 175, 200, 225 and 250 M non- 
rotating progenitors at Z = (the z-series) and Z = 
10~ 4 Zq (the u-series), which we summarize in Table [U 
The Z = 150 M star collapses to a black hole with- 
out exploding, so we exclude it from our CASTRO models. 
We employ metallicity rather than convective overshoot 
to cover the range of plausible progenitor structures in 
our study, because there is uncertainty about how much 
semi-convection there is in a given star but none about 
the range of metallicities over which it can form (Z = 
and 10~ 4 . Compariso n of Table □ with Table 1 of 
IScannapieco et aTl (j2QQ5h confirms that for a given pro- 
genitor mass these two metallicities do yield upper and 



lower limits to stellar radius similar to those for all rea- 
sonable values of semi-convection coefficients. 

3. RESULTS 

As expected, the u-series models die as red giants, with 
radii more than an order of magnitude larger than those 
of the z-series, which die as blue giants. As shown in Ta- 
ble [TJ the u-series models in general have more energetic 
explosions than the z-series models for a given mass. As 
we show in Figure [TJ convective mixing has completely 
disrupted the helium layer in model u225 and mixed it 
with the hydrogen envelope. The helium layer has also 
been mixed with the hydrogen envelope, although to a 
lesser extent, in model u200 and to an even slighter de- 
gree in model u250. These u-series models also experi- 
enced some mass loss due to pulsations. 

The pr 3 values through which the shock passes are an 
effective predictor of post-explosion hydrodynamics. In 
regions where pr 3 increases, the shock must decelerate, 
and a reverse shock forms that inverts the pressure gradi- 
ent and induces Ray leigh- Taylor (RT) instabilities. Mod- 
els with steeper values of pr 3 will experience more mixing 
because the forward shock slows down more abruptly, 
which leads to a stronger reverse shock. We show pr 3 
values for all models in Figure [TJ scaled to the maximum 
value in model u225. It is clear from this figure that 
models u225 and u200 exhibit the largest increases in 
pr 3 near the outer edge of the star. This is due to slight 
bumps in density that are connected to the pulsations 
that ejected mass from these stars earlier. 

We show final states of mixing for the PISNe in Figure 
[21 The models are shown at different times, but always 
after the shock emerged from the surface of the star and 
the reverse shock, if one formed, dissipated. Density is 
shown on the left in each panel, with values scaled to 
the maximum and minimum densities in the simulation, 
while oxygen is shown on the right, with the values in- 
dicated in the color bar. Some similarities between the 
models exist across all stellar structures. In the higher- 
mass compact models (z225 and z250) and in the red 
u-series models above 150 solar masses, the bulk of the 
oxygen layer was swept up into a shell of higher density 
than the material on either side. 

In general, however, the z-series explosions evolved 
quite differently than the u-series SNe. In the z-series, 
no reverse shock formed. A reverse shock formed in all 
the u-series models, but dissipated by the time it reached 
the density contrast at the oxygen layer (ul50 and ul75) 
or was too weak by the time it reached this dense shell 
to induce rapid RT growth (u250). In two models, a 
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Fig. 1. — The structure of helium and oxygen shells (delineated by dashed lines for clarity) in all the progenitors, with pr 3 superimposed 
on them. The pr 3 profiles have all been scaled to the maximum pr 3 in model u225. pr 3 increases dramatically near the edge of some 
models that die as red giants, while none of the blue giants show a similar structure. The plateaus in helium abundance above cosmological 
values, which are most apparent in model u225, denote convective regions in which the hydrogen envelope mixed with a portion of the 
helium layer. This process completely disrupted the helium layer in model u225, and mixed a smaller fraction of the helium layer with the 
hydrogen envelope in model u200. A small amount of the helium shell convectively mixed with the hydrogen envelope in model u250. 
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Fig. 2. — Images of density and oxygen abundance for all models after the forward shock exits the star or the reverse shock (if one forms) 
traverses the star, whichever is later. Density is scaled to the minimum and maximum values within a given simulation. A dense shell 
formed in all models with masses greater than 200 Mq, but the RT instability only grows appreciably in the u200 and u225 models. These 
were the largest progenitors in radius and had the steepest pr 3 curves. Only slight RT growth is evident in model u250. 
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reverse shock formed that was strong enough to drive 
rapid growth of RT instabilities (u200 and u225). The 
reverse shock was strongest in model u225, which had the 
steepest increase in pr 3 and the second highest explosion 
energy, after model u250. The RT instability is clearly 
visible in the u200 and u225 panels in Figure [21 In these 
models the instability became nonlinear because individ- 
ual RT fingers strongly interacted. In model u250, the 
RT instability only reached the early linear phase before 
its growth was halted. 

The extent to which the structure of the PISN at the 
time of explosion was preserved or disrupted by mixing 
is shown in Figure 02 We plot the initial shell structure 
of the PISN at the time of mapping into CASTRO as a 
dashed line, and overlay a solid line indicating the abun- 
dances of these elements after the shock exited the star 
or the reverse shock traversed the star, whichever was 
later. The slight smoothing in the final profiles in com- 
parison to the initial ones is due to numerical diffusion 
over the course of the simulation. The u200 and u225 
panels demonstrate the extent of mixing in these explo- 
sions. The oxygen layer has been completely disrupted, 
and mixing has penetrated just to the top of the silicon 
layer (in u200) or through the silicon layer (in u225). 
The most vigorous mixing occured in model u225, where 
hydrogen shell boost led to convection that completely 
mixed the helium shell and the hydrogen envelope. The 
only other models to experience RT mixing, models u200 
and u250, are also the only other models that show ev- 
idence of convective mixing in the outer region of the 
star prior to explosion. Model u200 manifests more RT 
mixing and a larger region (in radius and mass) in which 
convection mixed part of the helium layer with the hy- 
drogen envelope prior to explosion. Model u250, which 
has the smallest amount of convection extending into the 
helium layer, also exhibits the least RT growth. Models 
with no convective mixing between the helium layer and 
the hydrogen envelope manifest no RT mixing during the 
explosion, so mixing seems tied to the depth to which 
this convective envelope has penetrated the helium layer 
of the star. The deeper the convective envelope, the more 
RT mixing occurs during the supernova shock. 

RT-induced mixing in PISNe, if it occurs at all, is un- 
likely to leave the remnant in a state that resembles a 
core-collapse supernova remnant. In particular, it is un- 
likely to dredge up significant amounts of oxygen, let 
alone silicon or 56 Ni, into the outer regions of the star, or 
draw much hydrogen back towards the center of the blast. 
The most vigorous mixing occurs in the model with the 
helium shell that has been completely disrupted by con- 
vective mixing with the hydrogen envelope, and likely 
represents the most mixing that can occur in a PISN. 
Even in this model, 56 Ni cannot reach the upper layers 
of lighter elements. 

As noted earlier, our CASTRO simulations exclude the 
first stages of the explosion, in which oxygen burn- 
ing drives convect ive mixing that may perturb the star 
(|Chen et al.l l2Q10h . Since these perturbations are ex- 
pected to be of low mode and high amplitude, it is un- 
likely they would alter the essential conclusions of our 
survey. Reprising our two-dimensional calculations in 
three dimensions is unlikely to yield significantly differ- 
ent results. The RT instability initially grows about 30% 
faster in three dimensions than in two because of ar- 



tificia l drag forces arising from two-dimensional geome- 
tries ([Hammer et al.ll2010[ ). However, once the fingers of 
the instability begin to interact with one another, they 
mix more efficiently in three dimensions than in two, re- 
ducing the Atwood nu mber and hence their growth rate 
([Joggerst et al.l l2Q10bh . These two effects cancel each 
other, and the width of the mixed region in two and three 
dimensions is the same. For spherical simulations like 
ours in which the RT fingers significantly interact, two- 
dimensional and three-dimensional simulations would ex- 
hibit comparable degrees of mixing. Also, the RT insta- 
bility will not appear in three dimensions when it is not 
manifest in two. 

4. DISCUSSION AND CONCLUSIONS 

Unlike 15 - 40 M core-collapse Pop III SNe, 150 - 250 
M Pop III PISNe experience either no internal mixing 
prior to shock breakout from surface of the star or only 
modest mixing between the O and He shells. Minor mix- 
ing occurs in only two of the explosions and is due to 
the formation of a reverse shock that is strong enough to 
trigger the RT instability at the dense shell created by 
the forward shock at the top of the oxygen layer. The 
degree of mixing is principally a function of how well 
hydrogen shell boost mixes the helium shell with the hy- 
drogen envelope. The model in which the helium shell 
is completely mixed with the hydrogen envelope exhibits 
the most mixing. The general lack of internal mixing in 
PISNe has several consequences for early chemical en- 
richment of the IGM, the formation of second generation 
stars, and the observational signatures of such explosions. 

First, the elements that are expelled by low-mass Pop 
III SNe (which are governed by both mixing and fallback 
onto the compact remnant) are later imprinted on new 
stars in essentially the proportions in which they are cre- 
ated by the explosions, regardless of how and where the 
stars form. This is because the metals are already highly 
mixed by the time the shock exits the star and are merely 
diluted upon further expansion into the halo and IGM, 
where new stars might form. This implies that IMF aver- 
ages of nucleosynthetic yields of primordial core-collapse 
SNe can be directly compared to chemical abundances in 
ancient metal-poor stars, without regard for any inter- 
vening hydrodynamical processes. Indeed, the fact that 
elemental yields from Salpeter IMF averages of 15 - 40 
M progenitors in JET 10 match those found in two ob- 
servational surveys of extremely metal-poor (EMP) stars 
suggests that the bulk of early chemical enrichment may 
have been due to low-mass Pop III stars. This is at odds 
with the current state of the art in Pop III star forma- 
tion simulations, which suggest that the first stars were 
predominantly 100 - 500 M . 

Determining the nucleosynthetic imprint of PISNe on 
second-generation stars is much more problematic be- 
cause their metals may have differentially contaminated 
new stars. This is because metals in PISNe mixed with 
each other and with the surrounding IGM on much larger 
spatial scales, at radii where new stars may have formed. 
Except in models u200 and u225 at the interface of the 
O and He shells, and in model u225 at the interface of 
the O and Si shells, the shells of elements in the ejecta of 
other PISNe expand homologously until they sweep up 
their own mass in the ambient H II region, at radii of 10 
- 15 pc. At this point, a reverse shock detaches from the 
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Fig. 3. — Initial and final abundances of He, O, Si, and Fe for all models. The RT instability is only present in models u200, u225, 
and to a small degree in u250. The slight disagreement between initial and final profiles for models which experienced no mixing is due to 
numerical diffusion over the evolution times of these models. In none of the mixed models did mixing break through the Si layer of the 
star, so 56 Ni never reaches the surface of the explosion as it does in core-collapse supernova explosions. 
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forward shock and a contact discontinuity forms b etween 
them (the Chevalier phase- IWhalen et al.l 120081 ). Both 
the reverse shock and contact discontinuity are prone to 
dynamical instabilities that will mix elements from the 
interior of the remnant with the surrounding medium. 
Later, on scales of 100 - 200 pc, further mixing will oc- 
cur upon collision of the remnant with the dense shell of 
the relic H II region. 

In this case, the elements that are taken up into new 
stars depends on how and where the stars form. If met- 
als migrate out into cosmological flows and then fall back 
into the the halo via accretion and mergers on timescales 
of 50 - 100 Myr, they likely will be well-mixed, with all el- 
ements formed in the explosion appearing in the new sta r 
(jGreif et al .120071 IWise fe AbeHl2008bl : iGreif et al.ll20lol) . 
However, new stars may form at much earlier times in the 
SN remnant at radii where mixing takes place. One way 
this could happen is if metals at the interface between the 
ejecta and swept-up shell suffuse into and cool the shell, 
causing it to fragment into clumps that are un stable to 
gravitational collapse (e.g. iMackev et al.l 12003) . In this 
scenario, such clumps would be enriched only by the ele- 
ments residing in the relatively narrow zone in which the 
stars form. This still does not explain why some hyper 
metal-poor stars have such high [C, N, and O/Fe] ratios. 
While C and O would be preferentially deposited on the 
clumps in which such stars form because they are pre- 
dominant in the outer layers of the ejecta, PISNe do not 
produce enough N to account for measured abundances 
in these stars. However, little iron from deep in the in- 
terior of the remnant would reach the clumps because it 
is not mixed, which is consistent with the abundances 
measured in EMP stars to date. The failure to uncover 
the characteristic o dd-even nucleosyn t hetic signature of 
PISNe predicted bv lHeger fe Wooslevl (|2002l ) metal-poor 
star surveys thus far has led some to suppose that pri- 
mordial stars may not have been very massive, but it is 
possible that such signatu res have been masked by obser- 
vational selection effects (Ka rlsson et al. 2008). Unlike 
for low-mass stars, detailed numerical simulations that 
follow differential enrichment are required to determine 
the true chemical imprint of PISNe on second-generation 
stars. 

Second, our results imply that, unlike in SN 1987 A 
or low- mass Pop III SNe, Ni and Fe emission lines will 
not appear immediately after radiation breakout from 
the shock because there is not enough mixing to trans- 
port these elements out to the photosphere of the fireball. 
This may be key to distinguishing between Pop III core- 
collapse explosions and PISNe. The light curves of these 
supernovae are characterized by a sharp intense initial 
transient that decays over several hours into a dimmer 
extended plateau that persists for 2 - 3 mon ths in core- 
collapse events and 2 - 3 y ears in PISNe (|Frver et al.l 
120101 : IWhalen fe Frverl |2010D . The inital pulse is pow- 
ered by the thermal energy of the shock and the plateau 
is energized by radioactive decay of 56 Ni (the long life 
of the plateau is due to the longer radiation diffusion 
timescales through the massive ejecta). Preliminary ra- 



diation hydrodynamical calculations of PISN light curves 
and spectra indicate that their peak bolometric lumi- 
nosities are similar to those of Type la and core-collapse 
SNe, making determination of the mass of the progenitor 
from the magnitude of the initial transient problematic. 
If, however, Ni and Fe are detected just after radiation 
breakout, one can be confident that the explosion is due 
to a low-mass primordial star, and a Pop III IMF could 
begin to be built up by samples of such detections. 

Our models also demonstrate that one-dimensional ra- 
diation hydrodynamical models are sufficient to capture 
most features of PISN light curves and spectra because 
the mixing such calculations exclude, which would alter 
the order in which emission lines appear over time, is 
minor. The picture for core-collapse Pop III SNe is quite 
different because vigorous mixing prior to the eruption of 
the shock from the surface of the star mandates its inclu- 
sion in light curve models. Two-dimensional multigroup 
radiation hydrodynamical calculations of such spectra 
lie within the realm of current petascale platforms, but 
a less costly approach can incorporate mixing in one- 
dimensional models. If most mixing in low-mass Pop III 
explosions occurs before radiation breaks free from the 
shock, two-dimensional models such as those in JET10 
can be used to compute the distribution of elements in 
the ejecta just before breakout. These explosions can 
then be azimuthally averaged onto the one-dimensional 
grid of the light curve calculation and evolved to com- 
pute spectra. On average, along any given line of sight 
out of the SN, this method will produce emission lines in 
the likely order they would be observed. 

These simulations, together with our previous survey 
of mixing and fallback in low- mass Pop III SNe, are the 
first of a numerical campaign to model the chemical en- 
richment of the early cosmos from its smallest relevant 
spatial scales. The eventual goal of this campaign is to 
understand the contribution of the first SNe to the for- 
mation of new stars and the assembly of primeval galax- 
ies, which will soon be probed by the James Webb Space 
Telescope (JWST) and the Atacama Large Millimeter 
Array (ALMA). The next stage of these numerical sim- 
ulations will incorporate mixing on subparsec scales to 
determine if new stars directly form in the remnants of 
the first supernova explosions and follow their congrega- 
tion into the first galaxies. 
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